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ABSTRACT 

Chemical landscape of natural RNA species is 
decorated with the large number of modified nu- 
cleosides. Some of those could easily be detected 
by reverse transcription, while others permit only 
high-performance liquid chromatography or mass- 
spectrometry detection. Presence of m 6 A nucleo- 
side at a particular position of long RNA molecule 
is challenging to observe. Here we report an easy 
and high-throughput method for detection of m 6 A 
nucleosides in RNA based on high-resolution 
melting analysis. The method relies on the 
previous knowledge of the modified nucleoside 
position at a particular place of RNA and allows 
rapid screening for conditions or genes necessary 
for formation of that modification. 

INTRODUCTION 

More than a 100 modified nucleoside species could be 
found in various functional RNA molecules (1). 
Monitoring of their presence at a particular position of 
RNA is necessary for the identification of conditions af- 
fecting RNA modification as well as for the screening for 
genes responsible for the modification (2). Some of the 
modified nucleosides, e.g. m 2 G and m'A, impede reverse 
transcription and could easily be detected by primer ex- 
tension assay (3-6). Others need an additional chemical 
treatment of modified RNA before reverse transcription. 
For example, m 7 G nucleoside serves as a chain breakage 
point on treatment with NaBH 4 followed by incubation 
with aniline (7,8). Pseudouridine could be revealed by 
reverse transcription if RNA is pre-treated with N 2 H 4 or 



N-cyclohexyl-N'-(2-morpholmoethyl)carbodiimide 
metho-p-toluenesulfonate and is incubated in a mild basic 
solution (9). Adenosine nucleosides monomethylated at 
their exocyclic amino group (m 6 A) are among the 
hardest to detect. For monitoring of the m 6 A presence 
in a specified RNA position within a large RNA, a short 
fragment of that RNA should be excised. This could be 
done with either RNase H cleavage in the presence of the 
complementary oligodeoxyribonucleotide pair (10) or 
with unspecific RNase digestion of total RNA, while the 
selected region is protected by the complementary 
oligodeoxyribonucleotide (11). After purification of 
RNA fragment, it should be digested with the specific 
RNase, usually RNase A or RNase Tl. The resulting set 
of oligoribonucleotides should be analyzed by mass spec- 
trometry. Recently published alternative for this method is 
an immunoprecipitation with m 6 A-specific antibodies 
followed by a massively paralleled sequencing (12,13). 
Application of this method allowed revealing of thousands 
of m 6 A nucleosides in eukaryotic mRNAs. Monitoring of 
m 6 A presence at specific RNA positions in a variety of 
conditions is indispensable for understanding presumable 
regulatory role of m 6 A nucleosides. 

Here we report a new simple method based on high- 
resolution melting (HRM) analysis, which allows detec- 
tion of m 6 A nucleoside at the specific RNA position. 
The method could easily be adapted to a high-throughput 
format. 



MATERIALS AND METHODS 

Strains and media 

Escherichia coli strains BW25113 (WT) (14) and JW5107 
(AybiN), JW1860 (AyecP), JW4268 (AyjhP), JW0203 
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(AyafS), JW0904 (AsmtA), JW1123 (AymfD), JW4366 
(AlasT), JW2565 (AyfiF), JW2559 (AyfiQ, JW3900 
(AyiiV), JW1859 (AyecO), JW5501 (zl^g), JW0050 
(zlfog^), JW3466 (AyhiR), JW5543 (zljM/) and JW0200 
(AyafE) (15) were grown at 37°C in LB media to 0.6 A 60 o- 
Bacillus subtilis transformed with pKH80 plasmid coding 
for ErmC methyltransferase was grown similarly, with 
addition of 50 ug/ml ampicillin. 

HEK 293 cells {Homo sapiens) were cultivated at 37°C, 
5.5% C0 2 in Dulbecco's modified Eagle's medium:F12 
medium (Invitrogen) supplemented with 10% fetal 
bovine serum. 

Ribosome and RNA preparation 

Ribosomes and ribosomal subunits were prepared as 
described (16). Pure 23 S ribosomal RNA from the strains 
tested was prepared by phenol extraction out of the 50 S 
ribosomal subunits (17). Total RNA from bacteria was 
prepared by QIAGEN RNeasy mini kit. RNA species 
were stored in deionized water at — 80°C. To prepare 
total tRNA, SI 00 extract was made by the ultracentrifuga- 
tion of the cell lysate in MLA-130 rotor (80000rpm, 1 h). 



SI 00 was phenol extracted and then small RNA fraction 
containing all tRNA species was precipitated by 0.3 M 
NH 4 OAc and 80% ethanol solution. Each of RNA isola- 
tions was performed three times for three subsequent inde- 
pendent analyses (biological replicas). 

Total small RNA fraction (RNAs up to 200-nt long) 
from 10 8 HEK 293 cells was prepared by Ambion 
mirVana™ miRNA Isolation Kit. The approximate 
content of small nuclear RNA (snRNA) U6 was estimated 
by producing complementary DNA [cDNA; Maxima 
First Strand cDNA Synthesis Kit for Reverse transcript- 
ase-quantitative polymerase chain reaction (RT-qPCR), 
Thermo Scientific] and following RT qPCR, while the 
control during the amplification was U6 coding DNA 
with known concentration. 



Synthetic oligonucleotide probes carrying fluorophore and 
quencher, and snRNA U6 fragments 

We used following pairs of oligonucleotides (B-x 
quencher oligodeoxyribonucleotide, F-x-fluorescent oli- 
godeoxyribonucleotide): 



Name of oligodeoxy-ribonucleotide Sequence of oligodeoxyribonucleotide Complementary 

region of RNA 



Variation of quencher oligodeoxyribonucleotide length (A2030 of 23S rRNA E.coli). Figure 1A-D. 



B-2030 1 




5'-TCA CAG CG-3'-BHQl 


2023-2030 nt 


F-2030 1 




FAM-5'-AGT TCA ATT TCA CTG AGT CTC GG-3' 


2000-2022 nt 


B-2030 2 




5'-ATC TTC AC A GCG-3'-BHQl 


2023-2034 nt 


F-2030 1 




FAM-5'-AGT TCA ATT TCA CTG AGT CTC GG-3' 


2000-2022 nt 


B-2030 3 




5'-CTG CAT CTT CAC AGC G-3'-BHQl 


2023-2038 nt 


F-2030 1 




FAM-5'-AGT TCA ATT TCA CTG AGT CTC GG-3' 


2000-2022 nt 


B-2030 4 




5'-TACACTGCATCTTCACAGCG-3'-BHQ 1 


2023-2042 nt 


F-2030 1 




FAM-5'-AGT TCA ATT TCA CTG AGT CTC GG-3' 


2000-2022 nt 


Variation of fluorescent oligodeoxyribonucleotide length (A2030 of 23S rRNA E.coli). Figure 1E-H. 




B-2030 3 




5' -CTG CAT CTT CAC AGC G-3'-BHQl 


2023-2038 nt 


F-2030 5 




FAM-5'-AGT TCA ATT TCA CTG AGT-3' 


2005-2022 nt 


B-2030 3 




5'-CTG CAT CTT CAC AGC G-3'-BHQl 


2023-2038 nt 


F-2030 6 




FAM-5'-AGT TCA ATT TCA CTG AGT CT-3' 


2003-2022 nt 


B-2030 3 




5'-CTG CAT CTT CAC AGC G-3'-BHQl 


2023-2038 nt 


F-2030 1 




FAM-5'-AGT TCA ATT TCA CTG AGT CTC GG-3' 


2000-2022 nt 


B-2030 3 




5'-CTG CAT CTT CAC AGC G-3'-BHQl 


2023-2038 nt 


F-2030_7 




FAM-5'-AGT TCA ATT TCA CTG AGT CTC GGG T-3' 


1998-2022 nt 


Analysis of 


m 6 A methylation in 


A2030 of 23S rRNA E.coli. Figure 2A, 3-6. 




B-2030 8 




5'-TAC ACT GCA TCT T-3'-BHQl 


2030-2042 nt 


F-2030_8 




FAM-5'-CAC AGC GAG TTC AAT TTC ACT GA-3' 


2007-2029 nt 


Analysis of 


m 6 A methylation in 


A2030 of 23S rRNA E.coli. Figure 2B-D 




B-2030 9 




5'-GTA CAC TGC ATC TTC A-3'-BHQl 


2028-2043 nt 


F-2030 9 




FAM-5'-CAG CGA GTT CAA TTT CAC TGA GT-3' 


2005-2027 nt 


B-2030 2 




5'-ATC TTC AC A GCG-3'-BHQl 


2023-2034 nt 


F-2030 1 




FAM-5'-AGT TCA ATT TCA CTG AGT CTC GG-3' 


2000-2022 nt 


B-2030 10 




5'-CATCTTCACAGCGAGTTCAA-3'-BHQ 1 


2016-2035 nt 


F-2030_10 




FAM-5'-TTT CAC TGA GTC TCG GGT GGA GA-3' 


1993-2015 nt 


Analysis of 


m 6 A methylation in 


A37 of tRNA Vul , E.coli. Figure 7A-C 




B-37 1 




5'-GAC CCC CTC CTT-3'-BHQl 


37-48 nt 


F-37 1 




FAM-5'-GTA AGG GAG GTG CTC TCC CAG CT-3' 


14-36 nt 


B-37 2 




5'-CCC CCT CCT TGT A-3'-BHQl 


34-46 nt 


F-37 2 




FAM-5'-AGG GAG GTG CTC TCC CAG CTG AG-3' 


ll-33nt 


B-37 3 




5'-CTT GTA AGG-3'-BHQl 


32-39 nt 


F-37_3 




FAM-5'-GAG GTG CTC TCC CAG CTG AGC TA-3' 


8-31 nt 


Analysis of 


m 6 A methylation in 


A1618 of 23S rRNA E.coli. Figure 7D 




B-1618 




5'-GAC CAC CTG TGT-3'-BHQl 


1614-1625 nt 


F-1618 




FAM-5'-CG GTT TGG GGT ACG ATT TGA TGT-3' 


1591-1613 nt 


Analysis of 


m 6 A methylation in 


A2085 of 23S rRNA B.subtilis. Figure 7E 




B-2085 




5'-CGG GGT CTT TCC GT-3'-BHQl 


2081-2094 nt 



(continued) 
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Continued 



Name of oligodeoxy-ribonucleotide Sequence of oligodeoxyribonucleotide Complementary 

region of RNA 



F-2085 FAM-5'-CCT GTC GCG GGT AAC CTG CAT CT-3' 2058-2080 nt 
Analysis of m 6 A methylation in A43 of snRNA U6 H.sapiens. Figure 7F 

F-U6 FAM-5'-GTATCGTTCCAATTTTAGTATATGTGCTGC-3' 13-54 nt 
B-U6 5'-CTAATCTTCTCTG-3'-BHQl 
Fragments of U6 snRNA H.sapiens with and without m 6 A used for method calibration 

U6 + m 6 A 5'-GCUUCGGCAGCACAUAUACUAAAAUUGGAACGAUAC 7-65 nt 

m 6 dAGAGAAGAUUAGCAUGGCCCCUG-3' 

U6-m 6 A 5'-GCUUCGGCAGCACAUAUACUAAAAUUGGAACGAUAC 7-65 nt 

dAGAGAAGAUUAGCAUGGCCCCUG-3' 



Probes hybridization and HRM analysis 

Each hybridization and melting experiment was 
reproduced in three technical and three biological 
replicas (total nine replicas minimum). If not indicated 
otherwise, for hybridization and melting experiment, 
0.4 uM of individual RNA or RNA mixture was 
combined with 0.4-uM BHQ-x oligonucleotide and 
0.2 -uM FAM-x oligonucleotide in optically clear PCR 
tubes. Ionic conditions for the hybridization were 50 mM 
Tris-HCl, pH 8.3, and 40 mM KC1, the volume of mixture 
was 10 ul. The mixture was placed into the qPCR analyzer; 
in our case, we used CFX96 Real-Time System (Bio-Rad). 
The hybridization was done by heating to 80°C for 3 min 
followed by cooling to 20°C at a rate 3°/min. Melting curve 
were obtained by steady heating back to 80°C at a rate 1°/ 
min in increments of 0.5° and monitoring of FAM fluor- 
escence. Fluorescence curves were normalized and differ- 
ential curves were analyzed. The standard representations 
of differential fluorescence curves produced by HRM 
analysis software (18) were shown. 

To measure the suitability of melting analysis at various 
RNA concentrations, essentially the same protocol was 
used, whereas the concentrations of the 23S rRNA and 
oligonucleotides were taken as indicated in the text and 
figures. The total 23S rRNA concentration was kept 
constant at 0.4 uM, whereas the composition varied from 
100% of methylated 23S rRNA purified from the wild- 
type strain to 100% of unmethylated 23S rRNA purified 
from the JW3466 (AvhiR) strain (11) in increments of 
10%. For detection of tRNA Va \ modification at m 6 A37, 
4 tig of total tRNA was taken for each probe. For detec- 
tion of m 6 A presence in H. sapiens U6 snRNA, 1.4 tig of 
total small RNA fraction from HEK 293 cells was taken 
for each probe, which corresponds to 1.4pmol (0.14 uM) 
of U6 snRNA according to RT qPCR quantification. 
Same amount of pure synthetic analog of U6 snRNA 
fragment either containing m 6 dA43 or unmethylated 
dA43 was used as a control. We used m 6 dA instead of 
m 6 A because of lack of any commercially available m 6 A 
precursor for chemical synthesis of RNA. 

RESULTS 

Modified nucleosides in DNA and RNA could influence 
base-pairing or base-stacking interactions, which could 



alter melting properties of nucleic acid duplexes (19). We 
explored the possibility to use melting properties of RNA 
DNA duplex to monitor presence of modified nucleosides. 
We used the segment of E. coli 23S rRNA with the 
modified nucleoside m 6 A2030 as a first model. 
According to our previously published data, this nucleo- 
side is completely modified in the wild-type E. coli strain, 
while being completely unmodified in the JW3466 (AyhiR) 
strain, lacking the A2030-specific rRNA methyltransferase 
(11). To test the utility of HRM analysis for the monitor- 
ing of adenosine modification, we designed a set of 
oligodeoxyribonucleotide pairs complementary to the 
RNA target (Figures 1 and 2). One oligonucleotide in 
each pair was modified with 6-carboxyfluorescein (FAM) 
fluorophore at its 5'-end, whereas the other one carried 
Black Hole Quencher (BHQ1) quencher at its 3'-end. 
Sequences of both probes are designed in such a way 
that they should hybridize to target RNA side-by-side, 
bringing fluorophore and quencher to immediate juxta- 
position (Figures 1 and 2). Modified RNA nucleoside is 
always located within the borders of area covered with the 
shorter quencher-containing oligonucleotide. 

We tested how melting curves of BHQ1 -containing 
oligonucleotides are affected by A2030 nucleoside modifi- 
cation. These BHQ1 -containing oligonucleotides share the 
same 3'-end region and have variable length of 8, 12, 16 
and 20 nt. All of them are shorter than FAM-containing 
oligonucleotide, whose length was 23 nt (Figure 1). We 
assumed that fluorescence of FAM would be affected by 
the melting of the shorter BHQ1 -containing oligonucleo- 
tides, and this melting would be dependent on the presence 
of the modification in RNA. Obtained differential melting 
curves demonstrated dependence on the A2030 modifica- 
tion (Figure 1 B and C) in such a way that the modification 
lowered Tm of the shorter oligonucleotide. The recom- 
mendation regarding the length of BHQ1 -containing 
oligonucleotide is that it should be sufficiently large to 
effectively hybridize to the RNA, unlike 8mer oligonucleo- 
tide (Figure 1A), but it has to be as small as possible to get 
the maximal difference in Tm for the modified and un- 
modified RNA (see decrease in ATm in Figure IB to C 
to D). 

The length of the FAM-containing oligonucleotide was 
changed in the next set of the experiments (Figure 1E-H). 
We tested the 18-, 20-, 23- and 25-nt long FAM- 
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Figure 1. Optimization of probe length. Differential melting curves of duplexes formed by the 23S rRNA from the wild-type and AyhiR knockout 
strains of E. coli and a set of pairs of oligonucleotide probes. Schemes of oligonucleotide probes hybridization to the target-modified RNA molecule 
are depicted below each panel. Green curves correspond to the wild-type strain (23S rRNA modified at m 6 A2030), and red curves correspond to 
AyhiR knockout strain (23S rRNA unmodified at A2030). The length of probes varied as follows: (A) FAM probe 23 nt, BHQ1 probe 8nt; (B) FAM 
probe 23 nt, BHQ1 probe 12nt; (C) FAM probe 23 nt, BHQ1 probe 16nt; (D) FAM probe 23 nt, BHQ1 probe 20nt; (E) FAM probe 18 nt, BHQ1 
probe 16nt; (F) FAM probe 20nt, BHQ1 probe 16nt; (G) FAM probe 23 nt, BHQ1 probe 16nt; and (H) FAM probe 25 nt, BHQ1 probe 16 nt. 
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Figure 2. Optimization of probe location relative to modified nucleo- 
side. Differential melting curves of duplexes formed by the 23S rRNA 
from the wild-type and AyhiR knockout strains of E. coli and a set of 
pairs of oligonucleotide probes. Schemes of oligonucleotide probes hy- 
bridization to the target-modified RNA molecule are depicted below 
each panel. Green curves correspond to the wild-type strain (23S rRNA 
modified at m 6 A2030), and red curves correspond to AyhiR knockout 



containing oligonucleotides combined with the single 16- 
nt long BHQ1 -containing oligonucleotide. The shortest 
FAM-containing oligonucleotide demonstrated no de- 
pendence of Tm on the RNA modification (Figure IE). 
We assume that similar differential curves in this case 
result from the early melting of the duplex formed by 
the RNA and FAM-containing oligonucleotide, before 
the melting of the BHQ1 -containing oligonucleotide 
duplex, which is hybridized with the modification site of 
RNA. The longer is the FAM-containing oligonucleotide 
(see Figure IE to F to G to H), the greater is the difference 
in melting curves for modified and unmodified RNAs. 

To optimize hybridization position of BHQ1 -containing 
oligonucleotide relative to the site of the modification in 
RNA, we tested four pairs of BHQ1 and FAM-containing 
oligonucleotides (Figure 2). BHQ1 -containing oligo- 
nucleotides hybridized to different areas of RNA so that 
the modification site is either complementary to the 3'-end 
of the quencher oligonucleotide (Figure 2A), or to the 
middle part of the quencher oligonucleotide (Figure 2B 
and C) or to the 5'-end region of the quencher oligo- 
nucleotide (Figure 2D). We observed much more effect 
of the BHQ1 -oligonucleotide length rather than of the 
position of quencher oligonucleotide relative to the 
modified nucleoside in RNA. Curiously, the pair of oligo- 
nucleotides hybridized in such a way that modification site 
in RNA, located opposite to the site of BHQ1 attachment, 
produced the most pronounced difference between melting 
curves for modified and unmodified RNAs (Figure 2A). 
Duplexes formed by oligonucleotide probes with the 23S 
rRNA carrying methylated nucleoside A2030 melts via 
biphasic process with two melting peaks at 43° C and 
52° C, whereas the duplex formed by the unmethylated 
23S rRNA melts at 49°C (Figure 2A). Maximal difference 
between the methylated and unmethylated RNAs is mani- 
fested at 49°C. This pair of oligonucleotides (B-2030_8 
and F-2030_8) was chosen for several further experiments. 

Could the observation of melting curves of RNA with 
selected oligonucleotide pair be helpful in searching for 
gene responsible for particular RNA modification via 
screening a set of bacterial strains carrying gene knock- 
outs? To answer this question, we chose a set of strains 
each of which lacked particular gene coding for actual or 
hypothetical RNA methyltransferase, namely, AybiN, 
AyecP, AyjhP, AyafS, AsmtA, AymfD, AlasT , AyfiF, 
AyfiC, AyiiV, AyecO, AygiQ, AksgA, AyhiR, AyhdJ and 
AyafE strains (15). The wild-type strain (14) was used as 
a control. The 23S rRNA was extracted from these strains 
and was studied by HRM analysis for the m 6 A2030 
presence (Figure 3). In all the strains, except AyhiR 
strain (11), nucleoside A2030 is modified. Accordingly, 
we can see highly similar differential melting curves of 
the 23 S rRNA-oligo nucleotide duplexes for rRNA from 

Figure 2. Continued 

strain (23S rRNA unmodified at A2030). The length of probes varied as 
follows: (A) m 6 A nucleoside hybridize opposite to the 3'-end of BHQ1 
probe; (B) m 6 A nucleoside hybridize at a 3-nt distance from the 3'-end 
of BHQ1 probe; (C) m 6 A nucleoside hybridize at a 8-nt distance from 
the 3'-end of BHQ1 probe; (D) m 6 A nucleoside hybridize at a 15-nt 
distance from the 3'-end of BHQ1 probe. 
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Figure 3. Applicability of the method for screening a set of knockout strains in a search for specific RNA methyltransferase gene. Differential 
melting curves of duplexes formed by a pair of oligonucleotide probes, shown in Figure 2A, and the 23S rRNA from the wild-type and the set of 
E. coli knockout strains. Key to the figure is shown on the right. 



all of these strains, including wild-type rRNA. Only the 
duplex formed by the oligonucleotides and the 23S rRNA 
from the strain AyhiR, which lacks A2030 modification, 
displayed clearly altered melting properties (Figure 3). 

We checked if we could detect m A in modified RNA 
being mixed with unmodified RNA of the same sequence. 
To this end, we prepared a set of mixtures of the 23 S 
rRNA containing m 6 A2030 from the wild-type strain 
with that from AyhiR strain, devoid of the modification 
at A2030. We used a set of mixtures of modified and un- 
modified RNA in different proportions ranging from 
100% of modified RNA to 100% of unmodified in incre- 
ments of 10%. The total concentration of RNA was 
0.4 uM. The resulting HRM curves of the mixtures 
changed steadily from the 100% of methylated RNA to 
100% of unmethylated, displaying roughly an additive 
contribution of respective pure RNA species (Figure 4). 
Approximately 30% of admixture of RNA containing 
m 6 A2030 with unmodified one can be detected. 

Our next goal was to determine the minimal amount of 
the RNA, whose modification we could detect using 
HRM-based method. To do so, we tested the melting 
behavior of 23S rRNA:B-2030_8:F-2030_8 complexes 
with the concentration ratios of: 1 liM:1 liM:0.5 uM, 
0.4liM:0.4liM:0.2liM, 0.2 [iM:0.2 uM:0.1 uM, 0.1 yM: 
OA liM:0.05 liM and 0.05 uM:0.1 uM:0.025 uM. 0.02 llM: 
0.02liM:0.01 uM (Figure 5). The amplitude of the fluores- 
cence change decreased with the decrease of RNA concen- 
tration. The concentrations of 0.1 uM of 23S rRNA, 
0.1 (iM B-2030_8 oligonucleotide and 0.05 liM F-2030_8 
oligonucleotide were found to be minimal for the 
decisive identification of the RNA modification. In 10-^1 
volume, the quantity of RNA equals to 1 pmol, that 
corresponds to the amount of rRNA extracted from 
5 • 10 4 bacterial cells. The molar ratio of 0.4 liM:0.4 liM: 
0.2 liM seems to be optimal for clear implementing of the 
method. 



Hybridization of oligonucleotide probes to RNA is 
specific and thus the presence of bulk cellular RNA 
should not interfere with the detection of the modification 
at the specific site in the particular RNA. We compared the 
melting curve analysis done with the purified RNA from 
the large ribosomal subunit (Figure 6A) with that of total 
cellular RNA (Figure 6B). Use of total RNA proved to be 
as effective as usage of purified RNA in melting curve 
assay for the detection of the m 6 A2030 modification. 

To test the universality of the method, we used four 
additional models having the modified nucleoside in the 
different primary and secondary structure contexts. First 
is the E. coli tRNA Va \ that contains nucleoside m 6 A37 
modified by methyltransferase YfiC. We designed three 
pairs of oligonucleotides complementary to tRNA Val ! 
(Figure 7A-C) and investigated melting of their duplexes 
with total small RNA fraction from the wild-type and 
AyfiC strains (Figure 7A-C). Application of the 
quencher-containing B-37_l oligonucleotide that hybrid- 
izes with the region 37^8 nt of tRNA Va \ produced highly 
different melting curves for the RNA prepared from the 
wild-type and AyfiC strain. A distinct melting peak was 
observed for only unmodified tRNA Va \, whereas no 
melting was evident for the modified RNA sample 
(Figure 7A). This difference could be explained by 
overall lack of oligonucleotide hybridization to the 
modified tRNA Va \. Modified nucleosides in tRNA are 
known to stabilize its structure (20) and it is likely that 
such stabilization prevents efficient competition of com- 
plementary oligonucleotide with the secondary and 
tertiary structure of tRNA. Shifting oligonucleotide hy- 
bridization position by two nucleotides (Figure 7B) has 
not changed the principal difference between the curves 
corresponding to the modified and unmodified tRNAs. 
However, for both samples, a peak at similar higher tem- 
perature could be observed (Figure 7B, 72°C), which 
could represent either melting of a duplex formed by 
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FAM-containing oligonucleotide or other structural re- 
arrangement independent of the modification. Further 
movement of hybridization position resulted in complete 
loss of hybridization with both the modified and unmodi- 
fied tRNA (Figure 7C). Lack of hybridization could be 
explained by the smaller length of the BHQ-containing 
oligonucleotide. Alternatively, the presence of other 
modified nucleoside, uridine 5-oxyacetic acid, in the 
middle of the tRNA region complementary to the oligo- 
nucleotide prevents its hybridization. 

As a second additional model, we used the E. coli 23S 
rRNA carrying m 6 A modification at a nucleoside 1618, 
which is formed by YbiN methyltransferase. The pair of 
B-1618 and F-1618 oligonucleotides allowed to observe a 
reproducible difference between curves for the methylated 
23S rRNA from the wild-type and for the unmethylated 
23S rRNA from AybiN strain (Figure 7D). The difference 
observed for the m 6 A1618 modification resembled that 
for m 6 A2030 modification (Figure 1). In both cases, modi- 
fication resulted in lower melting temperature of the 
RNA-DNA hybrid. 

We make use of 23S rRNA adenosine 2058 (E. coli 
numbering) dimethylation introduced by ErmC 
methyltransferase, an enzyme that causes macrolide resist- 
ance in a number of bacteria (21) as our third additional 
model. Ribosomal RNA was purified from B. subtilis 
transformed with the plasmid coding for ermC gene. 
The similar RNA sample from untransformed B. subtilis 
was used as a control. We used pair of oligonucleotides B- 
2085 and F-2085 complementary to B. subtilis 23S rRNA 
(Figure 7E). The area of B-2085 oligonucleotide comple- 
mentarity in the 23S rRNA contained modified m 6 2 A2085, 
the nucleoside corresponding to E. coli m 6 2 A2058. 
Because adenosine dimethylation completely abolish 
Watson-Crick base pairing, we observed drastic difference 
between the melting curves corresponding to the modified 
and unmodified 23S rRNA (Figure 7E). Melting curve 
corresponding to unmodified 23S rRNA contain a single 



pronounced peak, while the melting behavior of the 
duplex formed by modified 23 S rRNA demonstrated 
several shallow peaks. 

Finally, we demonstrated that the method is suitable for 
the detection of m 6 A presence in eukaryotic non-ribosomal 
RNA. To this end, we used one of the few m 6 A residues in 
human RNA, which was previously identified at a single- 
nucleotide resolution, namely, m 6 A43 of U6 snRNA (22). 
Unfortunately, an enzyme responsible for this modification 
is unknown. Thus, we first compared model chemically 
synthesized fragments of U6 snRNA containing m 6 dA43 
or dA43. We used deoxynucleosides due to the unavailabil- 
ity of the m 6 A ribonucleotide precursor for the RNA 
chemical synthesis. After initial optimization of the fluor- 
escent and quencher oligonucleotide length, we chose the 
most optimal pair (Figure 7F). In complete agreement 
with our results on the previous models, we observed 
decrease in the melting temperature on m 6 A modification 
(Figure 7F, red curve-unmodified versus green curve- 
modified). Because the enzyme responsible for m 6 A43 for- 
mation in U6 snRNA is unknown, we compared a melting 
curve corresponding to the natural U6 snRNA present in a 
total small RNA fraction of human HEK 293 cells with 
that of control chemically synthesized U6 snRNA frag- 
ments. In agreement with our expectations, a melting 
curve corresponding to the natural human U6 snRNA 
(Figure 7F, blue curve) demonstrated even lower Tm 
than the methylated control RNA. The difference 
between the curves corresponding for the m 6 dA-containing 
model RNA and the natural U6 snRNA could be 
attributed to other modifications present in natural U6 
snRNA, such as several 2'OMe and ty. 



DISCUSSION 

Hereby we report a method for monitoring the presence of 
m 6 A nucleoside in a specific position of long RNA 
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Figure 5. Sensitivity of the method. Differential melting curves of 
duplexes formed by the 23S rRNA from the wild-type E. coli strain 
(solid curves) and AyhiR knockout strain (dashed curves) and a pair of 
oligonucleotide probes, shown in Figure 2A at different concentrations of 
the tested RNA and probes. The 23S rRNA:BHQ 1 -probe:FAM-probe 
ratios are as follows: 1 uM:l uM:0.5 uM (A), 0.4 uM:0.4 uM:0.2 uM (B), 
0.2uM:0.2uM:0.1uM (C), 0.1 uM:0.1 uM:0.05uM (D), 0.05 uM: 
0.1 uM:0.025 yM (E) and 0.02uM:0.02uM:0.01 uM (F). 



molecule using HRM analysis. Adenosine nucleosides 
methylated at their exocyclic amino group are hard to 
detect (3). Primer extension analysis could not be used 
for this purpose because m 6 A does not affect reverse tran- 
scription of RNA. RNA fragment excision and cleavage 
followed by mass-spectrometry is extremely laborious and 
could not be scaled up to high-throughput analysis. RNA 
immunoprecipitation via anti-m 6 A antibody followed by 
massively parallel sequencing allows transcriptome-wide 
identification of m 6 A nucleosides with 50-100 nt fidelity 
(12,13). However, it could be hardly used to screen many 
samples in a high-throughput assay. The method 
described here allows easy screening for the m 6 A 
presence at a particular position of RNA using total 
RNA sample and qPCR machine. Although the 
monomethylation of adenosine exocyclic amino group 
does not lead to alteration in hydrogen bonding within 
the nucleic acid duplex, it affects stacking interactions 
and accordingly its melting properties (19). 

Highly precise monitoring of fluorescence in the course 
of duplex melting is widely used for the monitoring poly- 
morphisms in DNA (23). It was also reported to be useful 
for studying RNA editing (24) and DNA methylation (25). 
Both methods described in the literature rely on the 
changes in nucleic acid sequence and as such are extensions 
of basic mutation detection technique based on HRM. 
Here we demonstrated the suitability of HRM method 
for the direct monitoring of the modified nucleoside 
presence in the RNA molecules. In contrast to the previ- 
ously published methods it does not rely on the change in 
nucleic acid sequence and does not involve any PCR-based 
amplification. The method described here relies directly on 
the alteration of nucleic acids duplex melting properties 
resulting from the RNA modification. 

According to our data, the presence ofm 6 A in the RNA 
part of the RNA-DNA duplex lowers its melting tempera- 
ture. The extent of Tm shift increases with the decrease in 
length of the quencher-containing oligonucleotide; 
however, the short quencher-containing oligonucleotide 
could fail to hybridize with RNA especially if the RNA 
target possesses strong secondary structure. We recom- 
mend 12-13-nt-long probes containing a quencher and 
>20-nt long probes, containing a fluorophore. We also 
could advise to design a quencher-containing oligonucleo- 
tide in such a way that after hybridization with RNA, the 
modified RNA nucleoside would be directly opposite to 
the 3'-terminal nucleotide carrying the quencher. 

We demonstrated the universality of the method. The 
detection of m 6 A nucleosides at three specific positions of 
rRNA, one position of tRN A and one position of snRNA 
was successful. Not only purified-specific RNA, but also a 
bulk cellular RNA is suitable for detection of m 6 A at a 
pre-defined position due to the specificity of oligonucleo- 
tide probe hybridization. Non-ribosomal targets need 
partial enrichment, such as simple molecular weight- 
based purification, easily accomplished by commercially 
available kits. 

Possible application of the method is the screening of 
the knockout/knockdown strain libraries in a search for 
the gene responsible for the formation of particular m 6 A 
nucleoside. Another possible application is the detection 
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Figure 6. Applicability of the method for bulk cellular RNA. Differential melting curves of duplexes formed by the RNA samples from the wild-type 
and AyhiR knockout strains of E. coli and a pair of oligonucleotide probes shown in Figure 2A. Green curves correspond to the wild-type strain (23S 
rRNA modified at m 6 A2030), and red curves correspond to AyhiR knockout strain (23S rRNA unmodified at A2030). (A) Purifieid rRNA from the 
large ribosomal subunit was used for hybridization with probes; (B) Bulk cellular RNA was used for hybridization with probes. 



of the particular m 6 A nucleoside presence at various 
growth or environmental conditions. The overwhelming 
majority of modified nucleosides in eukaryotic and par- 
ticularly human mRNAs are m 6 A nucleosides (12,13). 
Given particular types of mRNA are present at the 
quantity of up to 10000 per mammalian cell (26), 10 7 
cells would be required to analyze methylation of 
abundant mRNA with the help of HRM method. 
However, it is advisable to apply partial mRNA enrich- 
ment, such as affinity purification via polyA tails. The role 



of these modified nucleosides still remains mysterious, 
while some important regulatory role is likely. The 
proposed method might help to decipher dynamic events 
which lead to modification of particular RNAs. Before the 
application of the described method to eukaryotic mRNA, 
positions of the modified residues have to be mapped with 
the nucleotide resolution. We demonstrated the suitability 
of the method to study eukaryotic RNA species 
exemplified by U6 snRNA. The latter is the rare 
example of eukaryotic RNAs where precise position of 
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Figure 7. Applicability of the method for other m 6 A nucleosides in RNA. (A-C) Differential melting curves of duplexes formed by the tRNA Va 'i 
from the wild-type and AyfiC knockout strains of E. coli and a set of pairs of oligonucleotide probes. Schemes of oligonucleotide probes hybrid- 
ization to the target-modified RNA molecule are depicted below each panel. Green curves correspond to the wild-type strain (tRNA modified at 
m 6 A37), and red curves correspond to AyfiC knockout strain (tRNA unmodified at A37). (D) Differential melting curves of duplexes formed by the 
23S rRNA samples from the wild-type strain of E. coli and AybiN knockout strain. Scheme of oligonucleotide probes hybridization to the target- 
modified RNA molecule is depicted below the panel. Green curves correspond to the wild-type strain (23S rRNA modified at m 6 A1618), and red 
curves correspond to AybiN knockout strain (23S rRNA unmodified at A1618). (E) Differential melting curves of duplexes formed by the 23S rRNA 
samples from either untransformed B. subtilis or B. subtilis transformed with plasmid pKH80 encoding ErmC methyltransferase. Scheme of oligo- 
nucleotide probes hybridization to the target-modified RNA molecule is depicted below the panel. Green curves correspond to the untransformed 
strain (23S rRNA unmodified at A2085), and red curves correspond to the strain transformed by pKH80 (23S rRNA modified at m 6 2 A2085). (F) 
Differential melting curves of duplexes formed by the synthetic fragments of U6 snRNA containing either m 6 dA43 (green curve) or dA43 (red curve). 
Melting curve corresponding to natural human U6 snRNA purified from HEK 293 cell line is colored blue. Scheme of oligonucleotide probe 
hybridization to the target-modified RNA molecule is depicted below the panel. 



m 6 A has been mapped. Because an enzyme responsible for 
m 6 A formation in U6 snRNA is unknown, the designed 
system that could be used for the screening of putative 
methyltransferase knockout or knockdown cell lines in a 
search for such an enzyme. 

Another application of this method could be the moni- 
toring of modified nucleosides in rRNA that cause antibi- 
otic resistance. A number of m 6 A methyltransferases 
mono- or dimethylate the nucleoside A2058 of the 23S 
rRNA [for reviews see (27,28)]. This modification leads 
to the resistance to macrolide antibiotics (21). The 
method described in this work was tested for monitoring 
such modifications as well. Presence of other modified nu- 
cleosides whose presence in the RNA would cause alter- 
ation in melting could also be monitored by the method 
described. 
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